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ABSTRACT 


Measurements  were  made  at  three  levels  above  the  sea  sur- 
face to  observe  possible  wave  influence  on  turbulence  para- 
meters. Tvro  levels,  at  one  and  three  meters  above  the  sea 
surface,  were  mounted  on  a small  floating  platform  tethered 
to  the  bow  of  the  R/V  Acania.  The  third  level,  at  6.6  meters 
above  the  sea  surface,  was  mounted  on  a boom  on  the  bow  of 
the  ship.  Measurements  were  made  of  mean  wind,  mean  tempera- 
ture, mean  humidity,  and  of  the  fluctuating  temperatures  and 
velocities  at  all  three  levels.  These  data  were  interpreted 

on  the  basis  of  empirically  derived  expressions  describing 

2 

the  temperature  structure  function  parameter,  , and  the 

rate  of  dissipation  of  turbulent  kinetic  energy,  £.  For  opti- 

2 2 

cal  relevance,  CT  was  related  to  , the  refractive  index 
structure  function  parameter. 

Results  show  qualitatively  that  there  is  a substantial 
wave  influence  near  the  sea  surface,  and  that  the  system  and 
procedures  used  to  obtain  the  data  produce  valid  measurements 
for  deriving  the  wave  influence  on  turbulence  and  on  optical 
transmission  in  the  near  sea  surface  layer  adjacent  to  the 
waves . 
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I . INTRODUCTION 


Describing  near  sea  surface  variations  of  turbulence 
parameters  is  important  in  many  areas  such  as  the  coupling 
between  the  wind  and  the  waves  or  the  overall  heat  budget 
of  the  atmosphere  and  the  oceans.  In  this  study  the  primary 
concern  is  limited  to  the  effect  of  the  near  sea  surface 
variations  of  the  turbulence  parameters  on  photographic 
resolution.  In  particular,  the  effects  are  most  predominant 
when  the  line  of  sight  is  along  or  near  the  sea  surface,  as 
when  viewing  from  a small  craft,  or  when  the  line  of  sight 
is  over  the  horizon. 

Turbulence  parameters  investigated  in  this  study  were 

2 

Crp  , the  temperature  structure  function  parameter,  and  e , 

the  rate  of  dissipation  of  turbulent  kinetic  energy.  For 

optical  wavelengths,  , the  index  of  refraction  structure 

function  parameter,  has  been  shown  to  be  directly  related 
2 2 

to  CT  . Both  and  t are  related  to  image  resolution 
2 

since  describes  the  intensity  of  refractive  index  fluc- 

tuations and  e determines  the  micro-scale  of  turbulence. 
Previous  studies  by  Davidson,  et  al.  (1977)  of  over 

water  turbulence  parameters  have  established  the  height  varia- 
2 

tions  of  and  e for  levels  at  least  six  meters  above  the 

water.  However,  there  exist  few  data  to  verify  the  height 

2 

variations  of  C^,  and  e within  the  first  few  meters  above 
the  water.  It  is  the  layer  below  six  meters,  adjacent  to 
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the  sea  surface,  that  is  of  primary  interest  in  this  study. 

It  is  in  this  layer  that  the  wave  influence  on  turbulence 
parameters  is  expected  to  predominate. 

2 

Earlier  studies  have  shown  that  C^,  and  e vary  inversely 
with  height  depending  on  the  stability  of  the  atmosphere  in 
that  layer.  The  inverse  relationship  suggests  that  the 
higher  turbulence  can  be  expected  nearest  the  sea  surface. 

In  the  earlier  studies,  the  wave  influence  was  not  examined 
thoroughly.  It  is  suspected  that  the  wave  influence  on  the 
turbulence  adjacent  to  the  sea  surface  would  increase  the 
turbulence  over  waves  compared  to  the  turbulence  over  a flat 
sea  surface. 

The  nature  and  scope  of  this  study  is  primarily  to  estab- 
lish procedures  for  making  measurements  of  the  turbulence 
2 

parameters  and  e in  the  layer  adjacent  ot  the  sea  surface 

and  in  the  presence  of  some  wave  action,  mainly  swell.  Mea- 
surements were  made  from  about  six  meters  above  mean  sea 
level  to  within  one  meter  of  the  actual  sea  surface.  The 
levels  are  defined  in  reference  to  the  surface  of  the  longer 
period  swell  waves.  Secondarily,  the  results  of  these  mea- 
surements were  analyzed  to  see  if  any  wave  influence  could 
be  inferred,  and  if  so,  if  there  was  any  correlation  to  the 
variations  in  the  photographic  resolution  when  the  line  of 
sight  was  through  the  layer  adjacent  to  the  sea  surface. 
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II.  BACKGROUND 


A.  GENERAL  THEORY 

As  stated,  the  primary  concern  of  this  study  is  the  ef- 
fect of  the  variations  of  turbulence  parameters  on  the  photo- 
graphic resolution  when  the  transmission  path  is  near  the 
sea  surface.  Image  resolution  has  been  empirically  related 
to  the  scale  of  the  refractive  index  inhomogeneities  and  to 
the  intensity  of  fluctuations  in  the  refractive  index.  These 

in  turn  can  be  related  to  the  refractive  index  structure  func- 
2 

tion,  , which  can  be  directly  related  to  the  temperature 

2 

structure  function  parameter,  CT  , and  to  the  rate  of  dissi- 
pation of  turbulent  kinetic  energy,  e.  Similarity  expres- 

2 

sions  exist  which  relate  and  e to  mean  conditions.  The 

observation  objectives  of  this  study  were  to  measure  and 
analyze  the  mean  temperature,  wind  speed,  and  humidity  imme- 
diately adjacent  to  the  sea  surface. 

In  turbulence  theory,  it  is  important  to  consider  whether 
the  turbulence  is  isotropic  or  anisotropic.  In  general, 
large  scale  turbulence  is  anisotropic  and  small  scale  turbu- 
lence is  isotropic.  The  universal  formulae  considered  are 
valid  only  for  isotropic  small  scale  turbulence. 

Kolmogorov  (1941)  postulated  that  small  scale  turbulence 
should  be  isotropic  even  though  it  was  generated  by  and 
embedded  in  anisotropic  large  scale  turbulence.  If  the 
Reynolds  numbers  are  large  enough,  the  turbulence  will  adjust 
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through  inertial  transfer  and  viscous  dissipation  until  a 
statistical  equilibrium  is  reached. 

Since  the  small  scale  fluctuations  are  isotropic  in 
nature,  only  one  parameter  is  necessary  to  describe  the  in- 
tensity of  the  refractive  index  fluctuations  over  many 

scales.  It  is  the  refractive  index  structure  function  para- 
2 

meter,  , defined  as 

CN2  = [ <n(x)  - n(x+r)>]2/r2/3  (1) 

where  n(x)  and  n(x+r)  are  refractive  indices  at  two 

points  on  a line  oriented  normal  to  the  mean  wind  direction 

and  separated  by  the  distance  r , which  is  less  than  the 

outer  scale,  Lq  , the  lower  end  of  the  inertial  subrange,  and 

greater  than  the  inner  scale,  l , the  smallest  scale  of 

naturally  occurring  turbulence. 

The  refractive  index  is  determined  primarily  by  density 

fluctuations  and  can  therefore  be  related  to  temperature 

fluctuations.  Neglecting  humidity  effects  with  respect  to 

2 2 

optical  wavelengths,  the  expression  relating  to  is 
CN2  = [ 79. xlO-6  P/T2 ] 2 CT2  (2) 

The  expression  defining  the  temperature  structure  func- 

2 

tion  parameter  is  similar  to  the  one  for  and  is 

CT2  * [ <T ' (x)  - T’ (x+r)>]2/r2/3  (3) 
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where  T' (x)  and  T' (x+r)  are  temperature  fluctuations  at 

two  points  separated  by  the  distance  r . 

Practical  meteorological  and  optical  methods  for  direct 
2 2 

measurement  of  and  have  been  developed. 

2 

An  alternate  relationship  for  C^,  , based  on  measurements 
of  the  dissipation  of  turbulent  kinetic  energy,  e,  and  tem- 
perature variance,  x>  is 


(4) 


where  6 is  an  empirical  constant  with  a value  of  3.2.  This 

2 

form  enables  indirect  estimates  of  C^,  to  be  made  from  the 
mean  conditions,  since  z and  x are  easily  related  to 
boundary  layer  fluxes  and  to  profiles  when  steady  horizontal- 
ly homogeneous  conditions  exist. 

Small  scale  velocity  fluctuation  properties  are  also  of 
interest  in  optical  propagation  because  image  resolution 
has  been  empirically  related  to  the  inner  scale,  , which 
is  defined  as 


= (Y3/e)1/4  (5) 

where  y is  the  kinematic  molecular  viscosity.  e can  be 
obtained  from  either  one-dimensional  velocity  variance 
spectral  estimates  in  the  inertial  subrange  or  from  velocity 
structure  function  estimates. 
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Based  on  empirically  derived  expressions,  e and  x can 
be  functionally  related  to  mean  profile  and  flux  estimates 
(U*,  T*,  Z/L)  resulting  in  the  following  relationships 


and 


4>1  (Z/L) 


V*2 
X = VL 


4>2  (Z/L) 


(6) 

(7) 


L is  the  "Monin- Obukhov"  stability  length  and  is  defined 
as 


- 5 u*3 

L = — 

gk  9v'w 


(8) 


Businger,  et  al.  (1971),  empirically  determined  the  forms 
of  $^(Z/L)  and  <1> 2 (Z/L)  as 


(1  - 15  Z/L)1/4  - Z/L  , \ < 0 
1 + 3.7  Z/L  , | > ° 


0.74(1  - 9 Z/L) _1/ 2 
0.74  +4.7  Z/L 


(9) 


(10) 


where  -Z/L  and  +Z/L  correspond  to  unstable  and  stable 
conditions . 
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2 

It  is  desirable  to  have  expressions  which  relate  C^. 
and  at  some  level,  Z,  to  parameters  which  can  be  esti- 
mated from  more  readily  available  mean  parameters.  Such 
parameters  are  boundary  fluxes  (U*  and  T*)  or  mean  gradi- 
ents obtained  from  multiple  level  measurements  or  the 
wind  at  a given  level  and  an  air-water  temperature  differ- 
ence. Such  relationships  are 

CT2  = T*2Z"2/3  <j>3 (Z/L)  (11) 


where 


*3<r) 


*'2/3  ♦,<§ r1'3  *2<|) 


A formulation  for  <J>3(Z/L)  was  obtained  by  Wyngaard, 
et  al.  (1971)  on  the  basis  of  overland  data.  Their  formu- 
lation was 


$3(Z/L)  = 4. 9 [ 1-7 (Z/L) ] "2/3 

, 0 > Z/L 

(12) 

<t>3(Z/L)  = 4.  9 [ 1+2 . 75 (Z/L)  ] 

, 0 < Z/L 

(13) 

An  examination  of  Equations  (11) , (12) , and  (13)  with 

2 

the  knowledge  that  T*  is  independent  of  height  in  the  sur- 
face layer  yields  the  following  general  predictions  for  the 

2 

height  variations  of  CT  : 
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i)  For  near  neutral  conditions,  (Z/L  *0),  CT  will 

-2/3 

decrease  with  height  as  Z , since  Equation  (11) 


becomes 


m 2_2/3 


where  T*  is  independent  of  height. 

2 

For  unstable  conditions  (Z/L  <_  0)  , will  de- 


crease more  rapidly  with  height  than  Z 


-2/3 


-7  Z/L  >>  1,  C.p  should  decrease  with  height  as 


approximately  Z 
expressed  as 


-4/3 


since  Equation  (11)  can  be 


CT2  = 3 (££)2/3  (i5fT)4/3  z"4/3  (15) 


T 2/3  r 4/3 

where  (^)  (wT' ) is  independent  of  height. 

2 

iii)  For  stable  conditions  (Z/L  0)  , CT  will  decrease 

-2/3 

with  height  less  rapidly  than  Z . A limiting 
expression,  such  as  Equation  (13)  for  -7  Z/L  >>  1, 
is  not  considered  for  this  case  since  at  large 
positive  values  of  Z/L,  mechanical  turbulence  is 
nonexistent . 

The  inner  scale,  HQ,  can  also  be  related  to  the  boundary 
fluxes  and  mean  gradients  in  the  expression 


*Q  - [kZy3/U*3  ♦1(^)]3C 


r> 


An  examination  of  e on  the  basis  of  Equations  (6)  and 

3 

(9),  with  the  knowledge  that  U*  and  L are  independent  of 
height  in  the  surface  layer,  yields  the  following  general 
expressions  for  the  variation  of  e with  height. 

i)  For  neutral  conditions  (Z/L  = 0)  , e will  decrease 
with  height  as  Z ^ since  Equation  (6)  becomes 


•*  -1 
n 


ii)  For  stable  conditions  (Z/L  > 0)  , e will  decrease 
less  rapidly  with  height  than  Z”^  since  Equation 
(6)  becomes 


* -1 

Z La  + 3.  7 Z/L) 


and  (1+3.7  Z/L)  increases  with  height, 

iii)  For  unstable  conditions  (Z/L  < 0) , e will  also 

decrease  less  rapidly  with  height  than  Z~^,  since 
Equation  (6)  becomes 


e - -ir  £>“*- 


and  the  quantity  in  brackets,  (Z/L)  , increases 
with  increasing  height.  However,  this  increase 
with  height  is  less  than  for  the  stable  case,  so 
the  decrease  of  e with  height  in  the  unstable 
case  is  less  than  for  the  stable  case. 


1 


Selecting  the  mean  wind  at  two  appropriate  levels  (Z^ 
and  Z2)  and  when  conditions  are  at  or  near  neutral,  U*  can 
be  estimated  from  the  mean  wind  observation  as 

U*  - k(U2  - U1)/ln(Z2  - ZL)  (20) 

Furthermore,  when  conditions  are  at  or  near  neutral,  the 
production  of  turbulent  kinetic  energy  is  considered  to  be 
equal  to  the  dissipation  of  turbulent  kinetic  energy, 
yielding 

U*  * (ekZ)1/3  (21) 

The  relationships  given  by  Equations  (20)  and  (21)  allow  U* 
estimates  from  either  mean  wind  profiles  or  from  velocity 
fluctuation  data. 

2 

Alternate  expressions  for  and  e based  on  profile 

gradients  and  the  Richardson  number  (Ri)  as  the  stability 
parameter  are 


CT2  - z4/3oe/3Z)2  *4(Ri) 

(22) 

and 

e - (kZ)2O0/3Z)3J>5(Ri) 

(23) 

(g/T) ( 30/ 3Z) 

with 

Ri  - v 

(24) 

O0/3Z) 

where  $^(Ri)  and  ^(Ri)  are  empirically  derived  and  for 
mulated  on  the  basis  of  observations  of  temperature  and 
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velocity,  fluctuations  and  mean  gradients  of  wind  speed,  tem- 
perature, and  humidity,  and  where  9^  is  the  mean  virtual 
potential  temperature,  0 “ 8 (1  + 0 . 61  q)  . A formulation 

of  <t>^(Ri)  was  given  by  Wyngaard  et  al.  (1971). 

o 

A comparison  of  368  overwater  CT  and  profile  results 
with  the  overland  empirical  expression,  Equation  (22),  appears 
in  Figure  1.  The  distribution  predicted  by  Equation  (22) 
is  the  solid  line.  The  circles  are  mean  dimensionless  tem- 
perature structure  function  parameter  (DTSFP)  results  over 
.50  Richardon  number  (Ri)  intervals.  The  standard  deviation 
corresponding  to  each  mean  value  appears  as  a vertical  line 
on  which  the  number  of  values  defining  the  mean  appears. 

Larger  DTSFP  values  for  the  overwater  regime  for  near  neutral 
and  stable  conditions  (Ri  > -.5)  can  be  qualitatively  related 
to  the  wave  influence.  Decreased  DTSFP  values  for  unstable 
conditions  (Ri  < -.5)  are,  perhaps,  associated  with  aero- 


dynamic smooth  properties  for  the  sea  surface  during  light 
wind  conditions. 

2 

Evidence  of  wave  influence  on  C^.  appears  in  joint  proba- 
bility-conditional mean  results  of  the  normalized  temperature 
2 

variance  (a^  /T*)  over  waves  in  Figure  2.  The  results  and 
the  measurements  yielding  them  are  described  in  Davidson 
(1974).  The  ratio  (o^/T*)'^  is  directly  related  to  the  DTSP , 
Wyngaard,  et  al.  (1971).  In  Figure  2 the  vertical  axis 


variable  (C/U*)  defines  wave  influence  on  o.p/T*  , C is  the 
phase  speed  of  prevailing  swell,  L is  the  Monin-Obukhov  sta- 
bility length.  The  results  were  obtained  under  near  stable 
conditions.  The  increase  of  the  conditional  mean  of  o,p/T*  , 
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values  appear  as  numbers  on  the  grid,  for  increasing  C/U* 
for  near  neutral  conditions,  (Z/L  < .1),  provides  qualita- 


tive evidence  that  the  wave  influence  (C/U*)  on  o^/T*  , 
hence  the  DTSP , could  be  as  significant  as  the  stability 
influence  (Z/L) . 


%. 
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Figure  2.  Joint  probability-conditional  mean  results  of 
o_/TV  for  C/U*  and  Z/L  dependence.  (From 
Davidson,  1974) 

Volkov  (1969)  also  described  a wave  influence  on  o^/T*  . 
His  results,  Figure  3,  were  from  measurements  during  unstable 
conditions  and  established  a decrease  in  ctt/TVc.  due  to  the 
wave  influence,  C/U*.  This  results  agrees  with  those 
shown  in  Figure  1,  where  under  unstable  conditions  the 
DTSP  is  less  than  the  overland  prediction. 
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Figure  3.  Dependence  of  the  dimensionless  mean- square 
temperature  fluctuations  c_/T*  on  C/U*. 

(From  Volkov,  1969) 

A comparison  of  previously  observed  e results  and  the 
curve  based  on  Equation  (23)  appear  in  Figure  4 where  the 
predicted  curve  is  the  solid  line.  Averages  over  Richard- 
son number  (Ri)  intervals  of  .25  of  observed  data  appear 
as  small  circles.  The  error  bars  are  standard  deviations 
from  the  mean  within  each  interval  and  the  number  at  the  top 
of  each  error  bar  is  the  number  of  observations  defining  the 
mean  value.  The  significant  difference  between  the  predicted 
curve  and  the  observed  data  for  Ri  < .20  is  expected  based 
on  Eouation  (23).  For  Ri  > .20  the  equation  is  no  longer 
valid  since  terms  in  the  d>  ^ ( Ri ) expression  become  infinite. 

B.  PLATFORM  DESIGN 

Most  of  the  data  available  on  atmospheric  turbulence 
parameters  in  a marine  environment  were  obtained  from 


24 


sensors  mounted  on  ships  and  stationary  platforms  at  levels 
six  meters  and  higher  above  the  mean  sea  surface.  In  order 
to  detect  wave  influence  on  turbulence  parameters  and  to 
get  measurements  of  optical  transmission  in  the  layer  adja- 
cent to  the  waves  it  was  necessary  to  design  a system  with 
sensors  located  within  one  meter  of  the  actual  sea  surface. 

Several  designs  were  considered.  An  initial  design  was 
an  eight  meter  mast  mounted  on  a small  pontoon  type  float 
which  was  to  be  towed  astern  of  the  ship.  This  presented 
problems  relative  to  keeping  the  mast  nearly  vertical  yet 
have  the  float  small  enough  to  fit  on  the  after  deck  of 
the  R/V  Acania.  An  exceptionally  long  umbilical  cord  would 
also  have  been  required  so  the  platform  could  be  towed  far 
enough  astern  to  eliminate  ship's  influence  on  the  measure- 
ments. In  view  of  these  problems,  a system  was  selected 
which  could  be  deployed  from  the  bow  of  the  ship  and  could 
be  tethered  in  such  a manner  as  to  maintain  a position  in 
front  of  the  ship. 
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III.  DATA  COLLECTION 


A.  PLATFORM  AND  LOCATION 

All  observations  were  made  aboard  the  R/V  Acania  and 
from  a small  floating  platform  tethered  to  the  bow  of  the 
ship,  while  anchored  in  about  30  fathoms  of  water  off  Pt. 
Pinos  in  Monterey  Bay,  Figure  5.  Conditions  over  open 
ocean  differ  from  land  in  the  effects  of  wave  action  on 
turbulence,  in  the  nature  of  aerosols  and  fog,  and  humidity 
fluctuations.  These  conditions  could  best  be  obtained  far 
at  sea;  however,  the  cost  of  such  activities  made  it 
desirable  to  work  near  land  where  the  optical  propagation 
experiments  could  be  coincident.  Pt.  Pinos  and  other  loca- 
tions along  the  Monterey  Bay  shoreline  provided  very  nearly 
the  ideal  situations  for  both  the  optical  and  meteorological 
experiments . 

B.  FLOATING  PLATFORM  DESIGN 

The  final  design  utilized  a small,  light  weight  float, 

44  inches  wide  by  56  inches  long  by  12  inches  thick,  made 
of  styrofoam  and  plywood,  to  support  a ten  foot  hollow  mast, 
made  from  an  inch  and  one-quarter  mild  steel  pipe.  The  mast 
was  mounted  on  the  float  with  two  teflon  bearings  to  allow 
the  mast  to  turn  independently  of  the  float.  The  platform 
was  then  tethered  to  the  ship  by  passing  a steel  cable  down 
through  the  hollow  mast  and  out  through  the  bottom  of  the 
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Figure  5.  Optical  ranges  and  R/V  Acania  position  in 
Monterey  Bay. 


« 
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float.  A 100  pound  lead  weight  was  suspended  from  the  end 
of  the  steel  cable  to  give  enough  rigidity  to  the  cable  to. 
maintain  the  mast  in  a nearly  vertical  position.  The  steel 
cable  was  then  attached  to  a 15  foot  horizontal  bow  mounted 
boom  to  keep  the  platform  at  a position  ten  feet  in  front 
of  the  ship.  To  deploy  the  apparatus,  the  lead  weight  sus- 
pended from  the  tether  cable  was  lowered  to  at  least  ten  feet 
below  the  float  to  allow  for  vertical  movement  of  the  float 
as  it  followed  the  wave  surface.  The  float  was  designed 
with  enough  excess  buoyancy  so  that  it  would  follow  the  wave 
surface  with  very  little,  if  any,  lag  in  response. 

Two  sensor  levels  were  on  the  vertical  mast,  one  at  one 
meter  and  another  at  three  meters  above  the  water.  A wind 
vane  was  also  fixed  to  the  mast  to  keep  all  the  sensors 
pointed  into  the  wind.  A third  sensor  level  was  on  the  hori- 
zontal boom  above  the  float  approximately  seven  meters  above 
the  water.  The  umbilical  cord  ran  out  along  the  horizontal 
boom  and  then  to  the  top  of  the  mast.  A large  loop  main- 
tained in  the  cord  allowed  for  vertical  displacement  of  the 
float.  A small  hand  winch  mounted  on  the  horizontal  boom 
was  used  to  lift  the  entire  platform  with  the  steel  cable, 
so  it  could  be  swung  over  the  side  for  easy  deployment. 

Initial  dockside  tests  were  successful  and  only  minor 
difficulties  were  encountered  when  the  sensors  were  mounted. 
Initially,  a rudder  was  used  to  keep  the  float  from  spinning, 
but  this  interfered  with  the  sea-surface  temperature  probe. 

As  a final  modification,  the  rudder  was  removed  and  a small 
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wind  vane  was  mounted  on  the  float.  This  kept  the  float  from 
spinning  due  to  wave  action,  and  in  addition  reduced  the 
frictional  drag  on  the  mast  bearings  allowing  the  mast  wind 
vane  to  keep  the  sensors  more  closely  aligned  with  the  wind 
direction. 

Photographs  on  the  following  pages  show  the  platform, 
with  all  sensors  mounted  in  place,  being  deployed  from  the 
bow  of  the  R/V  Acania  with  the  ship  tied  alongside  the  Coast 
Guard  Pier  in  Monterey,  California.  There  was  no  wind  or 
wave  action  at  the  time.  The  platform  is  always  deployed 
while  the  ship  is  at  anchor.  However,  the  platform  could  be 
deployed  underway,  at  slow  speeds,  if  the  float  were  stream- 
lined, and  a streamlined,  heavier  lead  weight  were  used. 

C.  INSTRUMENTATION  AND  PROCESSING  EQUIPMENT 

Figure  8 shows  positions  of  the  instrumentation,  cable 
runs,  and  processing  and  recording  equipment  locations. 

The  mean  wind  measurements  were  obtained  with  cup  anemo- 
meter wind  profile  register  systems.  The  anemometer  consists 
of  three  cups,  plastic  cones  reinforced  with  aluminum  frames. 
When  the  cups  rotate  a slotted  shaft  in  the  anemometer  serves 
as  a shutter  between  a light  source  and  a photoelectric  cell. 
The  anemometers  have  a small  amount  of  internal  friction. 

Velocity  fluctuation  measurements  were  obtained  with 
hot  wire  anemometer  probes  with  sliding  support  shields  and 
tungsten  wires.  The  shields  permitted  isolation  of  the  sens- 
ing area  for  the  determination  of  the  undisturbed  velocity, 
VQ,  voltage  reading  before  and  after  each  experiment.  The 
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Figure  6.  Photograph  of  floating  platform  deployed  from 
the  bow  of  the  R/V  Acania. 
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platinum  coated  tungsten  wire  was  small  enough  to  resolve 
the  viscous  dissipation  scale  without  making  corrections 
for  wire  length.  The  anemometer  had  a linear  frequency  re- 
sponse from  0 KHz  (D.C.)  to  10  KHz  and  the  variable  decade 
module  operated  with  a 0-60  ohm  range. 

Temperature  fluctuations  were  measured  using  similar 
size  sensors.  A platinum  wire  was  used  in  place  of  a tung- 
sten wire.  The  temperature  fluctuation  system  was  designed 
for  a resolution  of  0.001°C  at  frequencies  up  to  1 KHz. 

The  temperature  measurements  require  a bridge  circuit. 

The  baseband  portion  of  this  bridge  circuit  is  basically 
a balanced  wheatstone  bridge  excited  by  a 3 KHz  signal  with 
a synchronous  detector  on  the  output.  The  very  small  dia- 
meter platinum  wire  serves  as  temperature  sensors  in 
opposite  arms  of  the  bridge.  The  resistance  temperature 
coefficients  result  in  an  output  from  the  bridge  which  is 
proportional  to  the  temperature  difference  between  the  two 

probes.  The  sensor  wire  is  0.5  centimeters  long  and  2 . 5x 
-4 

10  centimeters  in  diameter.  This  extremely  small  mass 
allows  a response  to  temperature  variations  of  up  to  1 KHz, 
while  electronic  amplification  allows  temperature  differ- 
ences as  small  as  0.004°C  to  be  observed. 

Both  wind  and  temperature  fluctuation  data  were  recorded 
on  a 14  channel  tape  recorder.  Real  time  readout  on  an 
eight  channel  chart  recorder  was  used  to  check  the  quality 
of  the  signals  coming  from  the  sensors.  The  chart  data  was 
also  utilized  in  temperature  variance  analyses. 
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Temperature  sensitive  quartz  crystal  probes  were  used 
to  obtain  mean  temperature  at  the  sea  surface  and  four  tower 
levels,  Figure  11.  The  RF  signal  from  the  probes  and  a 
reference  oscillator  were  mixed  in  a readout  unit  to  produce 
a beat  frequency  whose  signature  can  be  analyzed  to  within 
0.001  °C/Hz.  Each  sensor  received  pre-experiment  calibra- 
tion against  a platinum  resistance  wire  thermometer  in  a 
temperature  controlled  circulating  water  bath.  It  was 
checked  over  the  range  of  expected  temperatures.  The  accu- 
racy in  achieving  a 0.005°C  correction  factor  was  a constant 
for  each  probe.  The  tower  mounted  sensors  were  housed  in 
aspirated  shelters  to  eliminate  radiation  effects,  Figure  12. 

Mean  relative  humidity  information  was  obtained  using 
Dunmore  type  lithium  chloride  sensors,  Figure  13.  This 
sensor  was  also  placed  in  the  aspirated  shelter.  The  basic 
principle  of  operation  of  this  sensor  is  resistance  change 
in  an  electrolytic  solution  generating  a reference  voltage 
variation  proportional  to  relative  humidity  change.  Auto- 
matic temperature  compensation  in  the  instruments  meet  the 
following  specifications  for  relative  humidity,  + 3% 
humidity  below  90%  and  + 4%  humidity  above  90%.  Sensor  cali- 
bration was  accomplished  by  a comparative  method  using  a 
saturated  saline  solution  in  a closed  container. 

Sensor  placement  required  exceptionally  long  cable  runs, 
so  adjustments  were  made  in  the  bridges  for  resistance  and 
capacitance  of  the  wire  length  to  insure  good  frequency 
response . 
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thermometer  probe. 


Data  logging  during  the  experiments  was  accomplished 
using  the  NPS  developed  MIDAS  (Microprogrammable  Integrated 
Acquisition  System) . This  micro-processor  based  data  ac- 
quisition system  utilizes  an  Intel  8008  central  processor 
to  control  the  sampling,  averaging,  and  recording  of  mean 
meteorological  data.  All  software  is  written  in  PL/M  to 
facilitate  the  writing  of  a self-documenting  program. 

The  system  is  fully  automated  to  sample  the  tailored 
list  of  sensors  every  30  seconds  and  periodically  print 
output  values  averaged  over  the  selected  interval  of  from 
one  minute  to  one  hour.  In  this  study  a ten  minute  averag- 
ing interval  was  used.  Output  values  were  printed  out  on 
the  teletype  in  a columnized  format  with  the  time  of  print 
as  a leader.  A paper  punch  was  also  utilized  to  produce 
data  cards  for  use  with  the  IBM  360  digital  computer. 
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IV.  ANALYSES  PROCEDURES 

A.  ANALYSES  OF  MEAN  DATA 

Analyses  were  performed  on  data  recorded  on  teletype 
printout  which  was  controlled  by  MIDAS.  The  data  was  edited 
for  gross  errors  or  inconsistencies  due  to  known  instrumen- 
tation malfunctions.  The  criteria  for  retaining  or  discard- 
ing data  from  individual  levels  or  for  entire  ten  minute 
interval  depended  on  performance  check  results  obtained  dur- 
ing measurements,  obvious  inconsistencies  between  levels 
and  sequential  times,  and  the  availability  of  fluctuation 

I 

i data.  After  mean  data  were  obtained  for  ten  minute  inter- 

vals and  passed  preliminary  editing  for  obvious  erroneous 
values,  it  was  processed  for  three  applications,  i.e.,  mean 
wind  U,  mean  temperature  T,  and  mean  humidity  q.  These 
values  were  plotted  on  2-cycle  logarithmic  paper.  Since  U, 
T,  and  q are  expected  to  vary  logarithmically  with  height, 
a best  fit  straight  line  should  describe  the  distribution 
of  data  points.  Best  fit  lines  were  used  to  estimate  gra- 
dient parameters  and  the  Richardson  number. 

In  general,  the  procedures  were  subjective  and  in  some 
instances,  best  fit  lines  with  different  slopes  could  be 
visualized  for  a given  set  of  data  points.  Therefore,  one 
criterion  used  was  not  to  give  a single  data  point  too  much 
influence  in  determining  the  best  fit  line.  Consequently, 
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the  line  drawn  represented  a most  probable  position  between 
data  points  without  necessarily  coinciding  with  any  point. 

B.  ANALYSES  OF  PAIRED  TEMPERATURE  SENSOR 
2 

CT  values  were  estimated  by  using  paired  sensors  and 
variance  analyses.  This  analysis  was  based  on  the  expres- 
sion for  the  temperature  structure  function,  equation  (3). 

The  analysis  is  one  in  which  temperature  differences 
between  paired  resistance  wires  separated  by  a distance  r 
is  subjected  to  variance  analyses  as  indicated  by  the  term 
[ T ' (x)  - T'(x+r)]^  on  equation  (3). 

Voltages  corresponding  to  the  temperature  differences, 

AT,  as  measured  by  the  paired  sensors  were  recorded  on 
analog  magnetic  tape.  During  experimental  periods,  cali- 
bration voltages  were  also  recorded  on  the  data  tape  for 
reconstruction  of  the  signal.  This  calibration  signal  was 
played  back  at  the  time  of  data  analysis  to  provide  a 
reliable  check  of  the  recording  system  performance  in  both 
record  and  playback  modes. 

The  variance  or  mean  square  of  the  temperature  differ- 

2 

ences  required  to  compute  CT  was  obtained  by  analyzing  the 
analog  votages  with  a mean  square  voltmeter  using  a 30 
second  time  constant  for  the  averaging.  The  output  of  the 
mean  square  voltmeter  as  well  as  a record  of  the  original 
signal  were  recorded  on  standard  strip  charts.  A strip  chart 
section  appears  in  Figure  14.  The  record  of  the  original 
instantaneous  signal  provides  an  important  means  of  spotting 
extraneous  noises.  The  rms  voltmeter-strip  chart  recorder 
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system  was  also  calibrated  using  known  input  voltages.  At 
least  one  input  calibration  voltage  was  larger  than  the 
largest  expected  signal  on  the  data  record  and  likewise 
one  input  calibration  voltage  was  lower  than  the  smallest 
expected  signal.  This  procedure  minimized  any  non-lineari- 
ties in  either  the  rms  or  the  strip  chart  recorder. 

A final  check  on  the  validity  of  the  measured  data  was 
made  by  performing  spectral  analyses  on  a few  representa- 
tive data  samples  from  each  experimental  run.  The  spectrum 
analyzer,  described  in  section  C,  was  used  to  obtain  the 
spectra  which  were  subsequently  plotted  and  kept  as  a 
permanent  record.  Any  harmonic  noise  appeared  as  a spike 
in  an  otherwise  reasonably  smooth  curve.  It  should  be 
emphasized  that  continuous  calibration  and  data  checking 
procedure  were  essential  steps  throughout-  the  analyses. 

C.  ANALYSES  OF  VELOCITY  FLUCTUATION  DATA 

Spectral  analyses  to  obtain  e values  were  performed 
on  fluctuating  velocity  data  obtained  with  single  wires 
oriented  in  the  vertical  normal  to  the  mean  wind.  These 
data  represent  a point  measurement  and  yield  temporal  or 
time  descriptions  of  the  fluctuations.  Twenty-one  minute 
segments  of  data  originally  recorded  on  magnetic  tape  were 
recorded  into  a digital  spectrum  analyzer.  Procedures  for 
converting  spectral  values,  obtained  with  the  analog  spectral 
analyzer,  to  engineering  units  and  for  obtaining  turbulence 
parameters  from  the  spectra  are  described  in  the  following 
paragraphs . 
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A necessary  procedure  was  to  scale  the  spectral  plots 
in  order  to  relate  rms  input  voltages  to  power  spectral 
densities  (PSD) , variance  per  unit  frequency.  To  obtain 
PSD  levels,  corresponding  to  rms  voltage  input,  calibrated 
scale  charts  were  constructed. 

Amplitude  scaling  calibrations  were  accomplished  using 

externally  generated  white  noise.  Signals  with  a 1 volt  rms 

over  a frequency  range  of  0.1  Hz  to  1000  Hz,  giving  a PSD 
_ -3  2 

of  10  v /Hz.  Setting  a 3.16  v (10  dB)  input  gain  on  the 
spectrum  analyzer  insures  that  1 volt  input  corresponds  to 
Y=0.  An  example  of  such  a calibration  plot  is  shown  in 
Figure  15. 

e estimates  were  obtained  from  the  variance  spectra 
on  the  basis  of  the  universal  formulae  for  the  inertial 
subrange  in  wave  number  space.  These  expressions  predict 
a -5/3  slope  for  variance  spectra  of  both  velocity  and  tem- 
perature when  plotted  in  log- log  format.  Figures  16  and 
17  are  typical  spectra  of  both  variables  considered  in  the 
analyses.  The  velocity  variance  spectra  have  consistent 
-5/3  slope.  However,  temperature  spectra  often  exhibited 
slopes  slightly  different  from  the  expected  value  of  -5/3. 
This  feature  of  temperature  spectra  has  been  observed  by 
others  and  the  existence  of  "cold  spikes"  in  the  tempera- 
ture traces  has  been  postulated  as  a probable  cause  (Friehe, 
1976). 

Since  velocity  was  measured  at  a fixed  point  in  the  flow, 
the  resultant  spectral  estimates  correspond  to  "temporal" 
frequencies,  f.  To  obtain  e values,  temporal  (f)  and 
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Figure  15.  Calibration  plot  for  spectrum  analvzer. 
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Figure  16.  Example  of  temperature  variance  spectra. 
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Figure  17.  Example  of  velocity  variance  spectra. 
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space  (k)  scales  had  to  be  related.  This  was  accomplished 
by  Taylor's  (1938)  "frozen  turbulence"  hypothesis,  that  is 


k = 2Trf /U 


(25) 


where  U is  the  mean  wind  at  the  measurement  level.  The 
term  "frozen"  implies  the  turbulence  pattern  remains  un- 
changed as  it  sweeps  past  the  sensor  probe. 

The  following  identities  were  used  to  relate  temporal 
spectral  estimates  to  spatial  spectral  estimates 

f*Su(f)  = k-Su(k)  = Cl  e2/3  k_2/3  (26) 


where  S^(f)  is  the  spectral  density  value  with  units  of 
2 

(m/s)  /Hz.  The  above  equation  is  rewritten  in  the  follow- 
ing form  to  obtain  e from  values  of  f , U , and  ( f ) . 


e 


• f • S (f)}3/2 
u 


(27) 


where  k = 2irf/U  and  = . 5 and  C2  = .25,  empirically. 

D.  TEMPERATURE  SENSOR  AND  HOT  WIRE  CALIBRATION 

In-situ  calibrations  of  the  velocity  sensors  were  per- 
formed every  10  minutes  during  the  experiment.  For  this, 
recordings  were  made  of  both  the  cup  anemometer  wind  speed 
and  the  corresponding  hot  wire  voltage  output.  The  voltage 
from  the  sensor  is  given  by 

V2  - a 0*  + b (28) 
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where  V is  the  hot-wire  voltage  output  and  U is  the 
mean  wind  speed  for  any  given  level.  The  constants,  a and 
b,  are  the  calibration  curve  slope  and  intercept  respective- 
ly, and  were  derived  from  the  in-situ  calibration  curve. 

To  convert  the  voltages  to  velocity  units  requires  a 
calibration  factor  given  by 


U'  = Cv' 


(29) 


where  C is  the  calibration  factor  in  m/s  /volt.  v'  is 
the  voltage  fluctuation,  and  U'  the  velocity  fluctuation. 
Differentiating  equation  (28)  yields 


U'  = (4vU^/a)  v' 


(30) 


which  when  substituted  into  equation  (27)  yields 

C = 4vU*/a  (31) 

An  example  in-situ  calibration  curve  is  shown  in  Figure  15. 
The  slope  of  the  curve  yields  the  value  for  "a"  in  equation 
(31). 


51 


V.  RESULTS 


2 

Results  are  multi- level  mean  estimates  of  (derived 

2 

from  C^.  ) and  e . An  immediate  concern  in  interpretation 
was  the  possible  influence  of  the  ship's  structure  on  the 
measured  parameters.  The  absence  of  any  strong  ship  in- 
fluence can  be  established,  however.  In  Figure  18  are 
typical  profiles  of  the  mean  wind  data  obtained  during  this 
experiment.  The  mean  wind  versus  height  is  plotted  on  a 
log- linear  scale  producing  a straight  line,  which  is  to 
be  expected  in  the  neutral  case.  The  agreement  with  the 
neutral  case  prediction  for  the  mean  wind  profiles  suggests 
that  the  ship  influence  on  all  levels  was  negligible.  Other 
evidence  that  the  ship  influence  was  not  a factor  in  ob- 

2 

served  results  will  be  presented  in  the  discussion  of 
and  e results  which  follows. 

It  had  originally  been  planned  to  obtain  wave  spectra 
data  for  use  in  determining  the  swell  wave  speeds  in  the 
analyses.  However,  the  accelerometer  attached  to  mast  on 
the  floating  platform  did  not  function  properly,  and  no 
usable  interpretation  based  on  observed  wave  data  is  possible. 

(or  C.j.2)  Results 

The  primary  interest  in  this  study  was  to  examine  the 

2 

height  variation  of  the  optical  parameter  CN  . During  the 
entire  data  collection  period  the  stability  was  stable  or 
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very  near  neutral.  For  near  neutral  conditions,  neglect- 


ing the  wave  influence,  Equation  (11)  predicts  the  follow- 

2 

ing  height  dependence  for 


In  Figure  19  are  results  of  this  experiment  along  with 

the  results  from  earlier  experiments.  The  latter  are  shown 

in  order  to  emphasize  the  possible  slopes.  A best  fit 

line  to  the  upper  three  levels  closely  approximates  a -3/2 

slope,  agreeing  with  the  near  neutral  prediction. 

The  previous  examination  of  results  did  not  take  the 

wave  influence  into  account.  In  observational  studies  of 

the  wave  influence  on  turbulence  parameters,  Davidson  (1974) 

observed  the  wave  influence  to  be  as  significant  as  the 

stability  influence  as  the  sea  surface  is  approached.  This 

wave  influence  could  have  the  effect  of  increasing  the 

turbulence  in  the  air  immediately  above  the  waves  during 

2 2 

near  neutral  conditions  and  thus  increasing  CT  and 

above  the  waves  as  compared  to  over  a flat  sea  surface. 

A re-examination  of  the  results  in  Figure  19  at  the  6.6 

meter  level  and  below,  indicates  that  the  wave  influence  is 

2 

evident  by  a marked  increase  in  values  at  the  3 and  1 
meter  levels.  Even  though  conditions  were  near  neutral, 
the  profile  of  the  lower  three  points  approaches  the  rela- 
tionship for  the  unstable  case 


Qualitatively,  this  is  an  indication  that  the  wave  in- 
fluence very  near  the  sea  surface  may  override  the  stability 

2 2 

effect,  especially  on  CT  and  . 

Wave  influence  is  also  evident  in  the  near  neutral  por- 
tion (Ri  = 0)  of  the  summary  of  lower  level  results  data 
shown  in  Figure  20.  Here  it  can  be  seen  that  the  observed 
non-dimensional  overwater  results  in  general  have  higher 

values  than  those  predicted  based  on  only  stability. 

2 

In  addition  to  the  obtained  from  shipboard  meteoro- 

2 

logical,  CT  , data,  this  experiment  was  coordinated  with 

2 

a shorebased  laser  experiment  that  obtained  from  opti- 

cal data.  In  Figures  21,  22,  and  23  are  the  time  varia- 
2 

tions  of  at  the  three  levels  for  this  experiment  plotted 

2 

along  with  the  obtained  from  the  optical  experiment. 

The  results  labeled  MTF  are  from  optical  data  and  the 

results  BUS  are  from  meteorological  data  at  the  shore  based 

laser  site.  These  results  have  a much  closer  correlation 

2 

between  the  shore  b;.sed  optically  obtained  and  the 

2 

shipboard  meteorologically  obtained  than  previous  ex- 

periments have  shown.  This  closer  correlation  could  be 
due  to  the  fact  that  data  was  taken  at  this  low  level,  which 
was  very  near  the  mean  height  above  the  waves  of  the  laser 
beam's  transmission  path.  Previous  shipboard  measurements 
were  made  at  6.6  meters  and  higher, 
e and  U*  Results 

Figure  24  is  a plot  of  the  normalized  e results  versus 
height.  Again,  data  from  this  experiment  is  coupled  to  data 
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Figure  20.  Summary  of  results  of  low  level  DTSFP 
versus  Richardson  number. 
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2 

C„  versus  time  - shipboard  results  compared 
to  shore  results  for  2/15/77. 


12 


Figure  21. 
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Figure  23.  CM2  versus  tine  - shipboard  results  compared 
to  shore  results  for  2/17/77. 


from  earlier  experiments  to  emphasize  the  slope.  Except  for 
the  three  meter  level  value,  the  mean  profile  describes  the 
predicted  -1  slope. 

U*  results  for  this  experiment  are  shown  in  Figure  25. 
U*  is  plotted  for  each  level  that  dissipation  (e)  data 
was  available  and  is  also  shown  for  profile  (9U/9Z)  data. 

An  assumption  for  the  surface  layer  is  that  U*  is  constant 
with  height.  These  results  support  this  assumption,  in 
general,  since  U*  values  at  any  given  time  appear  to  be 
within  a factor  of  two  of  one  another.  Results  in  the  small 
table  in  the  upper  right  hand  corner  of  Figure  25  indicate 
that  the  mean  values  are  quite  close. 
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Figure  25.  Plots  of  U*  versus  height  for  each  level 
and  table  of  comparison  of  the  means. 


VI.  CONCLUSIONS 


The  purpose  of  this  study  was  to  use  previously  developed 
relationships,  based  on  recognized  sound  principles,  to 
validate  procedures  for  measuring  the  turbulence  parameters 
very  near  the  sea  surface.  The  experimental  results  do 
agree  reasonably  well  with  the  predicted  results,  substantiat- 
ing the  validity  of  the  new  procedures.  In  addition,  the 
results  showed  very  clear  evidence  that  the  wave  influence 
is  prominent  in  the  very  near  sea  surface  layer.  It  could  be 
as  important  a consideration  as  the  stability. 

There  is  more  development  needed  to  formulate  the 
theoretical  relationships  which  indicate  the  wave  influence 
on  the  near  sea  surface  turbulence  parameters.  This  experi- 
ment provided  a viable  and  reliable  method  for  collecting 
the  necessary  data  from  a ship.  The  use  of  the  tethered 
floating  platform  allows  expanded  use  of  already  proven  data 
collection  procedures  and  instrumentation,  thus  accomplish- 
ing the  primary  objectives  of  this  study. 
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APPENDIX  A 


DATA  TABLES 


TABLE  I 


DAY 

TIME 

TEMPERATURE  0 

T T T 

a 1 2 

G 

T3 

WIND  SPEED  M/SEC 

U1  u2  U3 

2/15/77 

H 

1540 

1600 

1620 

14.45 

15.14 

15.20 

15.25 

3.83 

— 

5.10 

• • 

14.35 

15.23 

15.37 

15.^9 

3.61 

— 

4.78 

• 1 

1640 

13.84 

15.10 

15.28 

15.44 

3.24 

— 

4.32 

II 

1700 

13.96 

14.84 

15.06 

15.25 

2.95 

— 

3.85 

II 

1720 

14.03 

14.56 

14.76 

14.93 

2.95 

— 

3.81 

2/16/77 

1350 

13.38 

12.93 

12.77 

12.59 

3.00 

— 

3.73 

II 

1400 

13.28 

12.82 

12.67 

12.52 

3.15 

3.51 

3.93 

ll 

1410 

13.25 

12.47 

12.34 

12.19 

3.34 

3.73 

4.22 

n 

1420 

13.23 

12.29 

12.15 

12.00 

3.39 

3.76 

4.24 

ii 

1430 

13.19 

12.21 

12.08 

11.94 

3.29 

3.71 

4.20 

ll 

1440 

13.16 

12.23 

12.12 

11.97 

3.76 

4.22 

6.76 

il 

1450 

13.13 

12.17 

12.07 

11.95 

4.32 

4.88 

5.61 

ll 

1500 

13.11 

12.12 

12.02 

11.91 

3.95 

4.42 

5.05 

ll 

1510 

13.08 

12.06 

11.96 

11.85 

3-78 

4.12 

^.73 

" 

1520 

13.07 

12.06 

11.97 

11.87 

3.29 

3.71 

4.05 

*• 

1540 

13.03 

12.03 

11.94 

11.83 

3.00 

3.2^ 

3.73 

•• 

1600 

13.00 

11.93 

11.85 

11.75 

2.73 

3.00 

3.34 

•I 

1620 

12.90 

11.96 

11.87 

11.77 

2.41 

2.68 

2.93 

•i 

1640 

12.95 

12.11 

12.04 

11.95 

2.51 

2.76 

3.07 

»• 

1700 

12.94 

12.29 

12.23 

12.15 

2.78 

3.07 

3.44 

* 

1720 

12.98 

12.45 

12.40 

12.3^ 

2.78 

3.39 

3.83 

1740 

12.98 

12.^9 

12.43 

12.35 

3.17 

. L'M. 

^.07 

** 

1800 

12.89 

12.53 

12.53 

12.52 

2.24 

2.46 

2.71 

2/17/77 

1420 

13.25 

14.83 

14.79 

14.73 

2.15 

2.46 

2.81 

II 

1440 

13.28 

14.82 

14.77 

14.71 

2.00 

2.27 

2.61 

II 

1500 

13.45 

14.88 

14.88 

14.86 

1.68 

1.98 

2.34 

1520 

13.45 

15.02 

15.02 

15.01 

1.68 

1.93 

2.29 

1540 

13.45 

15.33 

15.36 

15-38 

1.90 

2.12 

2.51 

1600 

13.51 

15.27 

15.23 

15. 2C 

1.93 

2.22 

2.51 

!• 

1720 

13.26 

14.65 

14.68 

1^.69 

1.63 

1.98 

2.> 

1740 

13.30 

14.14 

14.22 

14.29 

1.66 

1.98 

2.29 

•* 

1800 

13.36 

13.82 

13.92 

14. 00 

1.18 

1.95 

2.3^ 

1820 

13.^2 

13.72 

13.33 

13.93 

- • ' 

1.98 

2.36 

If 

1840 

13.40 

13.63 

13.74 

13.83 

1.73 

1.96 

2.37 

** 

1900 

13.41 

13.66 

13.75 

13.32 

1.95 

2.22 

2.71 

TABLE  II 


DAY 


2/15/?? 

H 


2/16/7? 

M 


2/17/7? 

M 


M 


TIME 

RELATIVE  HUMIDITY  % 
Hj_  H2  H3 

1540 

72.5 

73.3 

72.2 

1600 

— 

— 

— 

1620 

73.6 

71.2 

65.2 

1640 

73.2 

70.1 

64.8 

1700 

76.0 

71.7 

67.9 

1720 

79.2 

74.0 

72.1 

1350 

85.8 



89.6 

1400 

88.8 

86.0 

91.7 

1410 

91.2 

87.6 

93.4 

1420 

92.9 

88.7 

94.5 

1430 

94.1 

89.3 

95.1 

1440 

95.1 

90.1 

95.7 

1450 

95.4 

90.3 

96.2 

1500 

96.4 

91.4 

97.1 

1510 

97.4 

92.3 

97.7 

1520 

98.1 

92.9 

98.5 

1540 

98.6 

93,0 

99.0 

1600 

98.7 

93.0 

99.1 

1620 

99.5 

92.5 

99.7 

1640 

97.6 

93.0 

— 

1700 

— 

— 

— 

1720 

— 

— 

— 

1740 

— 

93.8 

94.0 

1800 

— 

93.5 

99.1 

1420 

71.8 

72.3 

72.4 

1440 

75.3 

74.0 

72.7 

1500 

75.8 

72.8 

70.3 

1520 

73.7 

71.2 

66.1 

1540 

70.5 

67.5 

63.6 

1600 

72.3 

70.0 

66.0 

1620 

— 



— 

1720 

78.3 

78.1 

74 .? 

1740 

85.7 

84.1 

81.0 

1800 

90.6 

88.8 

85.0 

1820 

92.3 

89.7 

85-7 

1840 

92.6 

90.2 

86.3 

1900 

92.9 

90.5 

86.6 

SPECIFIC  HUMIDITY  gm/kg 


*0 

ql 

q2 

q3 

10.31 

7.81 

7.93 

7.84 

10.24 

7.98 

7.79 

7.22 

9.91 

7.86 

7.63 

7.12 

9.99 

8.04 

7.69 

7.37 

10.04 

8.22 

7.87 

7.67 

9.62 

8.02 



8.19 

9.56 

8.28 

7.90 

8.34 

9.54 

8.27 

7.8  7 

8.31 

9.53 

8.32 

7.8? 

8.36 

9.50 

8.40 

7.89 

8.32 

9.48 

8.49 

7.98 

8.40 

9.4? 

8.48 

7.97 

8.43 

9.45 

8.54 

8.04 

8.48 

9.43 

8.60 

8.09 

9.50 

9.43 

8.66 

8.15 

8.58 

9.40 

8.69 

8.14 

3.60 

9.39 

3.70 

8.05 

9.57 

9.33 

8.73 

8.06 

3.63 

9.36 

8.64 

8.19 

9.37 



8.49 

8.46 

9.32 

— 

8.51 

9.02 

9.54 

7.59 

7.63 

7.61 

9.56 

7.96 

7.79 

7.63 

9-67 

8.04 

7.72 

7.45 

9.67 

7.89 

7.62 

7.07 

9.6? 

7.70 

7.39 

6.97 

9.71 

7.86 

7.60 

7.15 

9.55 

8.18 

8.20 

7.82 

9.58 

8.17 

8.55 

3.27 

9.61 

8.17 

8,85 

3.52 

9.65 

9,08 

5.89 

8.55 

9.63 

9.06 

8. 89 

8.55 

9.64 

9.11 

5.92 

3.58 
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TABLE  III 


DAY 

TIME 

CN 

S 

CO 

1 

O 

-2/3 

e 10  ^ m^/sec^ 

30/3Z 

10  u* 

1 

2 

3 

1 

2 

3 

10  c 

m/sec 

2/15/77 

1540 

25. 

9.3 

6.6 

14.0 

2.4 

1.1 

+0.87 

.190 

1600 

24. 

8.4 

5.5 

10.0 

1.3 

1.5 

— 

1620 

25. 

7.9 

5.1 

13.0 

1.1 

1.0 

+1.57 

.180 

1640 

26. 

7.9 

5.5 

5.6 

1.0 

0.92 

+2.06 

.160 

1700 

24. 

7.1 

4.3 

5.5 

1.1 

0.74 

+2.45 

.140 

1720 

20. 

5.7 

2.6 

5.0 

0.81 

O.63 

+2.25 

.130 

2/16/77 

1350 

21. 

7.4 

2.9 

5.3 

0.49 

0.43 

-1.50 

.110 

1400 

21. 

7.0 

2.9 

6.2 

0.49 

0.38 

-1.30 

.120 

1410 

20. 

7.0 

4.0 

6.3 

1.0 

0.57 

-1.19 

.130 

1420 

20. 

7.0 

4.3 

6.8 

1.7 

O.83 

-1.24 

.130 

1430 

20. 

6.6 

2.3 

6.2 

0.90 

0.76 

-1.13 

.140 

1440 

20. 

7.0 

2.3 

13.0 

1.7 

1.1 

-1.08 

.150 

1450 

20. 

6.1 

2.4 

13.8 

1.1 

0.96 

-0.87 

.190 

1500 

18. 

7.0 

2.1 

10.0 

0.99 

0.80 

-0.82 

.160 

1510 

19. 

6.1 

2.1 

8.0 

0.95 

0.78 

-0.82 

,140 

1520 

17. 

5.5 

1.7 

6.2 

0.90 

0.75 

-0.71 

.110 

1540 

19. 

4.8 

— 

4.8 

0.85 

0.42 

-0.77 

,110 

1600 

20. 

4.1 

— 

5.0 

O.63 

0.67 

-0.66 

,090 

1620 

19. 

3.7 

— 

— 

— 

— 

-0.71 

.055 

1640 

17. 

4.0 

1.4 

5.2 

0.72 

0.64 

-0.55 

,083 

1700 

13. 

2.2 

1.7 

4.0 

0.62 

0.66 

-0.45 

.098 

1720 

11. 

1.7 

2.2 

6.5 

0.90 

0.57 

-0.29 

.160 

17^0 

8.5 

1.8 

2.3 

4.3 

0.72 

0.81 

-0.44 

.130 

1800 

8.1 

1.3 

3.1 

3.7 

0.39 

0.15 

-0.24 

.070 

2/17/77 

1420 

3.2 

3.7 

1.8 

2.0 

0.45 

0.58 

-0.24 

.099 

1440 

3.5 

3.7 

1.1 

1.8 

0.30 

0.20 

-0.29 

.091 

1500 

3.8 

4.1 

— 

1.7 

0.31 

0.10 

-0.19 

.098 

1520 

— 

4.8 

2.3 

1.2 

0.04 

0.17 

-0.23 

.090 

1540 

— 

4.1 

2.8 

1.7 

0.064 

0.35 

+0. 56 

.090 

1600 

— 

4.1 

2.0 

4.0 

0.40 

0.46 

+0.50 

.037 

1620 

— 

— 

— 

— 

— 

— 

— 

1720 

5.2 

1.5 

0.77 

2.3 

0.28 

O.O32 

+0.50 

.098 

1740 

4.4 

1.4 

1.1 

2.5 

0.15 

0.014 

+1.08 

.094 

1800 

5.2 

1.6 

1.0 

2.6 

0.28 

0.014 

+1.24 

,098 

1820 

3.5 

1.8 

1.0 

2.6 

0.28 

0.027 

+1.40 

,088 

1840 

4.0 

1.9 

1.1 

2.8 

0.30 

0.010 

+1.34 

,095 

190C 

3.c 

1.9 

1.1 

4.6 

0.43 

-1.14 

,110 
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TABLE  IV 


DAY 

TTJ-E 

**10 

TlO 

°C 

U10 

m/sec 

9* 

2/15/77 

1540 

+.10 

15.25 

5.33 

— 

1 1 

1600 





— 

It 

1620 

+.12 

15.52 

4.99 

.38 

1 1 

1640 

+.21 

15.49 

4.51 

.39 

1 f 

1700 

+.39 

15.32 

4.02 

.34 

It 

1720 

+.42 

15.00 

3.96 

.27 

1 1 

1740 

— 

— 

— 

2/16/77  1350 

-.51 

12.49 

3.89 

_ 

1 1 

1400 

-.39 

12.43 

4.07 

— 

f t 

1410 

-.29 

12.11 

4.37 

— 

1 1 

1420 

-.30 

11.92 

4.39 

— 

l f 

1430 

-.24 

11.86 

4.37 

— 

1 1 

1440 

-.20 

11.98 

4.93 

.04 

It 

1450 

-.10 

11.88 

5.84 

.04 

it 

1500 

-.13 

11.84 

5.24 

.02 

ll 

1510 

-.17 

11.78 

4.90 

.02 

tl 

1520 

-.24 

11.81 

4.18 

.05 

II 

1540 

-.26 

11.76 

3.86 

.04 

II 

1600 

-.34 

11.96 

3.45 

.04 

II 

1620 

-.97 

11.71 

3.00 

.06 

II 

1640 

-.32 

11.89 

3.17 

.05 

ll 

1700 

-.19 

12.10 

3.56 

— 

tt 

1720 

-.05 

12.29 

4.02 

— 

tl 

1740 

-.11 

12.30 

4.22 

— 

II 

1800 

-.20 

12.50 

2.79 

— 

2/17/77  1420 

-.10 

14.69 

2.93 

_ _ _ 

II 

1440 

-.30 

14.66 

2.72 

.16 

II 

1500 

-.05 

14.83 

2.46 

.30 

ll 

1520 

-.16 

14.99 

2.40 

.40 

1 1 

1540 

-.06 

15.37 

2.62 

.36 

1 1 

1600 

-.40 

15.19 

2.61 

.38 

1 1 

1620 

— 

— 

— 

II 

1720 

+.07 

14.68 

2.46 

.18 

1 1 

1740 

+.33 

14.30 

2.40 

.20 

II 

1800 

+.36 

14.03 

2.46 

.22 

tl 

1820 

+.48 

13.96 

2.47 

.27 

II 

1840 

+.40 

13.85 

2.48 

.25 

II 

1900 

+.22 

13.83 

2.84 

.26 
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TABLE  V 


Profile  Level  1 Z=1.0m. 


DAY 

TIME 

AT 

X10 

ae/azio 

u* 

Ri 

u* 

r 

G 

io~2 

10"2 

2/15/77 

1540 

-.80 

-.28 

+ .87 

.19 

+.010 

.17 

8.60 

10.0 

ll 

1600 

— 

— 

— 

— 

.15 

— ~ 



II 

1620 

-1.37 

-.55 

+I.67 

.18 

+.012 

.17 

2.30 

3.50 

It 

1640 

-1.65 

-.81 

+2.06 

.16 

+ .021 

.13 

1.60 

2.60 

It 

1700 

-1.36 

-.84 

+2.45 

.14 

+.039 

.12 

1.00 

3.20 

II 

1720 

-0.97 

-.62 

+2.65 

.13 

+ .042 

.12 

0.82 

4.30 

2/16/77 

1350 

+O.89 

+ .60 

-1.50 

.11 

-.051 

.13 

2.00 

5.80 

• l 

1400 

+O.85 

+ .51 

-1.30 

.12 

-.037 

.13 

2.70 

6.30 

ll 

1410 

+1-.14 

+ .60 

-1.19 

.13 

-.029 

.13 

2.90 

3.20 

ll 

1420 

+1.31 

+.  68 

-1.24 

.13 

-.030 

.13 

2.70 

2.40 

ll 

1430 

+1.33 

+ .70 

-1.13 

.14 

-.024 

.13 

3.30 

2.40 

ll 

1440 

+1.18 

+.49 

-1.08 

.15 

-.023 

.17 

3.60 

3.00 

ll 

1450 

+1.25 

+.37 

-0.87 

.19 

-.010 

.17 

5.50 

2.70 

ll 

1500 

+1.27 

+ .46 

-0.82 

.16 

-.013 

.15 

5.00 

2.10 

ll 

1510 

+1.30 

t.54 

-0.82 

.14 

-.017 

.14 

5.60 

2.20 

ll 

1520 

+1.26 

+ .72 

-0.71 

.11 

-.024 

.13 

6.00 

1.90 

it 

1540 

+1.27 

+.8  5 

-0.77 

.11 

-.026 

.12 

6.30 

2.30 

ll 

1600 

+1.31 

+1.10 

-0,66 

.09 

-.034 

.12 

9.60 

2.40 

ll 

1620 

+1.13 

+1.26 

-0.71 

.055 

-.097 

7.50 

2.90 

II 

1640 

+1.06 

+1.06 

-0.55 

.083 

-.032 

.12 

10.0 

2.70 

ll 

1700 

+0.84 

+0.66 

-0.4  5 

.098 

-.019 

.11 

8.70 

2.50 

ll 

1720 

+O.69 

-•-0.43 

-0.29 

.16 

-.005 

.13 

15.0 

2.60 

il 

1740 

+0,68 

+O.38 

-0.44 

.13 

-.011 

.11 

3-90 

1.60 

ll 

1800 

+0.39 

+0.50 

+0.24 

.070 

-.020 

.11 

12.0 

4.50 

2/17/77 

1420 

-1.44 

-1.68 

-0.24 

.099 

-.010 

.089 

1.10 



H 

1440 

-1.38 

-1.86 

-0.29 

.091 

-.030 

.086 

— 

— 

• l 

1500 

-1.38 

-2.28 

+0.19 

.093 

-.005 

.084 

— 

— 

ll 

1520 

-1.54 

-2.67 

+0.23 

.090 

-.016 

.075 

— 

— 

ll 

1540 

-1.86 

-2.71 

tO.  56 

.090 

-.006 

.084 

— 

— 

ll 

ll 

1600 

-1.37 

-2.01 

+0.50 

.087 

-.040 

.11 

— 

— — 

II 

1720 

-1.42 

-2.35 

+0.50 

.098 

+ .007 

.093 

1.10 

1.40 

II 

1740 

-1.00 

-1.74 

+1.08 

.094 

+.033 

.096 

0.17 

0.20 

II 

1800 

-0.84 

-1.06 

+1.24 

.098 

+.036 

.097 

0.18 

0.69 

ll 

1820 

-0.54 

-0.91 

+1.40 

.083 

+.048 

.097 

O.O65 

0.44 

It 

1840 

-0.45 

-0.73 

+1.34 

.095 

+.040 

.099 

0.092 

0.82 

ll 

1900 

-0.17 

-0.21 

+1.14 

.11 

+ .072 

.12 

0.072 

3.20 

where  G - C ?2  Z4/3/AT2 

f - Ct2/Z4/3(36/3Z)2  - DTSFP 
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